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Kazuo Hida*, Ken'ichi Takano 1 , and Hidenori Suzuki 1 ^ 



o 

(N 

> 

O 

o 



in 

C3 



Division of Material Science, Graduate School of Science and Engineering, 

Saitama University, Saitama, Saitama, 338-8570 

1 Toyota Technological Institute, Tenpaku-ku, Nagoya 4-68-851 

(Received November 11, 2010) 

The ground states of two types of distorted mixed diamond chains with spins 1 and 1/2 
are investigated using exact diagonalization, DMRG, and mapping onto low-energy effective 
models. In the undistorted case, the ground state consists of an array of independent spin- 
1 clusters separated by singlet dimers. The lattice distortion induces an effective interaction 
between cluster spins. When this effective interaction is antiferromagnetic, several Haldane 
phases appear with or without spontaneous translational symmetry breakdown (STSB). The 
transition between the Haldane phase without STSB and that with (n + l)-fold STSB (n = 1, 
2, and 3) belongs to the same universality class as the (n + l)-clock model. In contrast, when 
the effective interaction is ferromagnetic, the quantized and partial ferrimagnetic phases appear 
with or without STSB. An effective low-energy theory for the partial ferrimagnetic phase is 
presented. 

KEYWORDS: mixed diamond chain, distortion, frustration, Haldane phase, spontaneous translational sym- 
metry breakdown, partial ferrimagnetism 
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1. Introduction 

Quantum magnetism in frustrated spin systems is a 
rapidly developing field of condensed matter physics. 1, 2 > 
At first glance, one would expect that geometrical frus- 
tration enhances quantum fluctuation and drives an or- 
dered state into a disordered state. However, recent 
progress in this field of physics has shown that this sim- 
ple intuition is not always valid and that geometrical 
frustration induces a variety of exotic quantum phenom- 
ena, which arc not easily predicted. Under an appropriate 
condition, it even stabilizes an unexpected magnetic long 
range order such as the frustration-induced ferrimagnetic 
and spin nematic orders. 

To understand magnetism under the interplay of geo- 
metrical frustration and quantum fluctuation, it is desir- 
able to begin with typical spin models with exact solu- 
tions. Among them, there exist a class of models whose 
ground states are exactly written down as spin cluster 
solid (SCS) states because of frustration. A SCS state is 
a tensor product state of exact local eigenstatcs of cluster 
spins. Well-known examples are the Majumdar- Ghosh 
model 3 ' 1 whose ground state is a prototype of sponta- 
neously dimerized phases in one-dimensional frustrated 
magnets 4 ) and the Shastry-Sutherland model 5 ' 1 which 
corresponds to the material SrCu2 (1303)2 - 6,7 ' 1 In these 
models, the spin clusters arc singlet dimers. 

The diamond chain is another frustrated spin chain 
with exact SCS ground states. The lattice structure is 
shown in Fig. 1. fn a unit cell, there are two kinds of 
nonequivalcnt lattice sites occupied by spins with mag- 
nitudes 5 and t; we denote the set of magnitudes by (5, 
t). One of the authors and coworkers 8, 9 ) introduced this 
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Fig. 1. Structure of the diamond chain. Spin magnitudes in a unit 
cell are indicated by S and r; we denote the set of magnitudes 
by (S, t). The PDC is the case of S = T, while the MDC is the 
case of S = 2t with an integer or half-odd integer r. 



lattice structure and generally investigated the case of 
(5, 5), i.e., the pure diamond chain (PDC). Any PDC 
is shown to have at least one exact SCS ground-state 
phase where each spin cluster has spin 0. Particularly, 
in the case of (1/2, 1/2), they determined the full phase 
diagram of the ground state by combining rigorous ar- 
guments with numerical calculations. After that, Nigge- 
mann et al. 10, n > argued about a series of diamond chains 
with (5, 1/2). As for the special case of (1/2, 1/2), they 
reproduced the results of ref. 9. 

The mixed diamond chain (MDC) is defined as a dia- 
mond chain with (5, 5/2) for the integer 5. 12 ) The spe- 
cial case of (1, 1/2) was first investigated by Niggemann 
et al. 10, n ) They considered it as one of the series of mod- 
els with (5, 1/2). Recently, extensive investigation on the 
MDC has been carried out by the present authors. 12 ~ 14 ) 
The MDC is of special interest among diamond chains, 
because only the MDC has the Haldane phase in the ab- 
sence of frustration, so that we can observe the transition 
from the Haldane phase to a SCS phase induced by frus- 
tration. In contrast, diamond chains of other types have 
ferrimagnetic ground states for weak frustration. 

The features common to all types of diamond chains 
arc their infinite number of local conservation laws and 
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more than two different types of exact SCS ground states 
that are realized depending on the strength of frustra- 
tion. For example, S = 1/2 PDC has a nonmagnetic 
phase accompanied by spontaneous translational symme- 
try breakdown (STSB) and a paramagnetic phase with- 
out STSB. This model also has a ferrimagnetic ground 
state in the less frustrated region. 9 ' On the other hand, 
the MDC with spins 1 and 1/2 has 3 different param- 
agnetic phases accompanied by STSB and one paramag- 
netic phase without STSB. This model also has a non- 
magnetic Haldane ground state in a less frustrated re- 
gion. 10,12 ) The SCS structures of the ground states are 
also reflected in characteristic thermal properties, as re- 
ported in ref. 13. 

Modifications of the PDC and MDC have been exam- 
ined by many authors. Among them, the spin 1/2 PDC 
with distortion has been thoroughly investigated by nu- 
merical methods. 15 ~ 17 ' It is found that azurite, a nat- 
ural mineral, consists of distorted PDCs with spin 1/2 
and that the magnetic properties of this material have 
been experimentally studied in detail. 18,19 ' Other mate- 
rials have also been reported. 20,21 ' The diamond chain 
is one of the simplest models compatible with the 4-spin 
cyclic interaction. The effects of this type of interaction 
on PDC have recently been investigated by Ivanov et 
al. 22 ^ The present authors also investigated the MDC 
with bond-alternating distortion and found an infinite 
series of ground states with STSB. 14 ' In addition, as re- 
viewed in ref. 14, the MDC is related to other impor- 
tant models of frustrated magnetism such as the dimer- 
plaquctte model, 23 ~ 28 ' frustrated Heisenberg ladders, 29 ' 
hybrid diamond chains consisting of Heisenberg bonds 
and Ising bonds, 30,31 ' and an Ising model on a hierarchi- 
cal diamond lattice. 32 ' Among them, the dimcr-plaquettc 
chain with ferromagnetic intcrplaquctte interaction re- 
duces to the MDC in the limit of strong interplaquettc 
interaction. 28 ' 

Thus far, in spite of the theoretical relevance of the 
MDC, no materials described by the MDC have been 
found. Nevertheless, synthesizing MDC materials is not 
an unrealistic expectation in view of the success of the 
synthesis of many low dimensional bimetallic magnetic 
compounds 33 ' and organic magnetic compounds. 34 ' In 
general, it is natural to expect that the lattice is pos- 
sibly distorted in real MDC compounds as in azurite. 
From this viewpoint, it is important to present theoret- 
ical predictions on the ground state of distorted MDCs 
to widen the range of candidate materials of MDC and 
to raise the possibility of their synthesis. 

We begin by classifying the distortion patterns by 
the normal modes of each diamond unit. Excluding two 
translations and one rigid body rotation, we have 5 nor- 
mal modes as depicted in Fig. 2 within the diamond 
plane. A distorted MDC may be realized as a result of 
the collective softening of these normal modes. In partic- 
ular, the distortion patterns in (a) and (b) break the lo- 
cal conservation laws that hold in the undistortcd MDC. 
Hence, these distortions induce effective interactions be- 
tween the cluster spins in the whole lattice, and may 
form novel exotic phases. We investigate these interest- 
ing cases in the present paper. In what follows, we name 
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Fig. 2. Displacement modes of a diamond unit. 
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Fig. 3. Structures of MDC with S = 1 and rW = r( 2 > = 1/2 
with (a) type A and (b) type B distortions. 



the distortion patterns in (a) and (b) as type A and type 
B, respectively. The MDCs with type A and type B dis- 
tortions are depicted in Figs. 3(a) and 3(b), respectively 
The distortion patterns in Figs. 2(d) and 2(e) do not 
change the geometry of the original undistorted MDC. 
The distortion pattern in Fig. 2(c) is of another interest, 
since it induces the bond alternation in the undistortcd 
MDC without breaking the local conservation laws. This 
case has been investigated separately and published in a 



previous paper 



14) 



This paper is organized as follows. In §2, the Hamil- 
tonians for the MDCs with the type A and type B dis- 
tortions are presented, and the structure of the ground 
states of the MDC without distortion is summarized. The 
ground-state phases for the MDC with the type A distor- 
tion are discussed in §3, and those for the MDC with the 
type B distortion are discussed in §4. The last section is 
devoted to summary and discussion. 

2. Hamiltonian 

The MDCs with the type A and type B distortions are 
described, respectively, by the following Hamiltonians: 

N 

H A =Y / [(1 + Sa)S 1 t{ 1) + (1 - 5 a )t, (1) 5,+i 



i=i 



(1 - S a )S 1 t\ 2) + (1 + <5a)t, (2) ^ +1 + AT, (1) T, (a) 



(1) 
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N 

Hb=J2 [(1 + Sb)S 1 t\ 1) + (1 + &)t, (1) Si+i 



1=1 



_(2) 



(1 - S b )Sit\" + (1 - fe)r, w S, +1 + Ar ; (1) x ; (2) 



.(2) 



(2) 



(2) 



where Si is the spin-1 operator, and t\ and t\ are 
the spin- 1/2 operators in the Zth unit cell. The param- 
eter 6a (5b) represents the strength of type A (type B) 
distortion, and is taken to be nonnegativc without spoil- 
ing generality. The number of unit cells is denoted by 
N, and then the total number of sites is 3N. We will 
consider these systems in the large N limit. 

For <5a = and 5b = 0, both eqs. (1) and (2) reduce 
to the undistorted MDC Hamiltonian, 



"Ho 



N 

E 

i=i 



SiTi + TiSi+i 



3 

2 



(3) 



.(2) 



with the composite spin operators T; 
Before going into the analysis of the distorted MDC, 
we briefly summarize the ground-state properties of the 
Hamiltonian (3) reported in ref. 12 for convenience. 

(i) T; commutes with the Hamiltonian Ho for any I. 
Therefore, the composite spin magnitude 7] defined 
by Tf = Ti(Ti + 1) is a good quantum number that 
takes the values or 1. Hence, each energy eigen- 
state has a definite set of {Ti}, i.e. a sequence of 
0's and l's with length N. A pair of r, and t\ 
with Ti — is called a dimer. A cluster including 
n successive 7] = 1 pairs bounded by two T; = 
pairs is called a clustcr-n. The clustcr-n is equiva- 
lent to an antiferromagnetic spin-1 Heiscnbcrg chain 
of length 2n+l with open boundary condition. Since 
a cluster-n is decomposed into a sublattice consist- 
ing of ti + 1 sites with Si 's and that consisting of n 
sites with T/'s, the ground states of a cluster-n are 
spin triplet states with total spin unity on the basis 
of the Lieb-Mattis theorem. 35, 36 ) This implies that 
each cluster-n carries a spin-1 in its ground state. 

(ii) There appear 5 distinct ground-state phases called 
dimcr-cluster-n (DCn) phases with n = 0,1,2,3, 
and co. The DCn state is an alternating array 
of dimers and cluster-n's. The phase boundary 
X c (n,n') between DCn and DCn' phases are 

A c (0,l) = 3, 

A c (1,2) -2.660425045542, 

A c (2, 3) -2.58274585704, 

A c (3,oo)~ 2.5773403291, (4) 

where A c (0, 1) is obtained analytically and other val- 
ues are calculated numerically. 

(iii) In the DCoo ground state realized for A < A c (3, oo), 
Ti = 1 for all I. This state is not accompanied by 
STSB and is equivalent to the Haldane state of an 
antiferromagnetic spin-1 Hcisenberg chain with infi- 
nite length. 

(iv) Each of the DCn states with < n < 3 realized for 



A > A c (3,oo) is a uniform array of cluster-n's with 
a common value of n and dimers in between. In the 
DCn phase with 1 < n < 3, (n + l)-fold STSB takes 
place. In the DC0 phase, no translational symmetry 
is broken. 

In what follows, we numerically examine various aspects 
of the type A and type B distortion effects on the MDC. 
Because the DC3 phase is only realized within a very 
narrow interval of A, it is difficult to analyze the effect 
of distortion numerically in this phase. Hence, we do not 
consider the DC3 phase in the following numerical anal- 
ysis. 

3. Ground-State Properties of the MDC with 
Type A Distortion 

3.1 Weak distortion regime 

We now inspect the nature of the effective interaction 
between two cluster-n's separated by a dimer consisting 



of 



.(i) 



and r, in the presence of the weak type A dis- 



tortion. For 5\ > 0, Si (Si+i) tends to be antiparallcl to 



.(i) <„m 



(2) /,.(!)> 



r, (t ) rather than to r, (t\ '), as is known from 

Fig. 3(a). The spins r, and r, are antiparallel to each 
other because they form a singlet dimer. Therefore, Si 
and Si+i tend to be antiparallel to each other. In each 
cluster-n, the number of spins Si's is larger than the 
number of composite spins TVs by one. Hence, from the 
Lieb-Mattis theorem, 35 ) the total spin of the ground state 
of the cluster-n points to the same direction as the Si 's 
belonging to that cluster-n. Therefore, the total spins of 
the cluster-n's on both sides of the dimer also tend to be 
antiparallcl to each other. Thus, the effective coupling 
between the spins of neighboring cluster-n's is antifer- 
romagnetic. This physical argument will be numerically 
ensured below. 

In general, the interaction between two spins with a 
magnitude of 1 is the sum of bilinear and biquadratic 
terms. Therefore, the effective Hamiltonian for cluster- 
n's in the phase that continues to the DCn phase in the 
limit of 5a — > is written as 



N c 



n cS = ^H off (i,i + i), 



(5) 



i=l 



H en {i,i + l) = J e &S i S i+ i+K eS {S i S i+1 ) , (6) 

where Si is the total spin of the i-th cluster-n with a 
magnitude of 1, N c is the total number of cluster-n's, 
and J e ff and K e s are effective coupling constants. From 
symmetry consideration, the signs of 5a does not affect 
the sign of the effective coupling constants. Hence, these 
coupling constants are of the order of 5 A for small 5a- 
We numerically calculated the ground-state energy of a 
pair of cluster-n's with total spin Stot, an d compared 
it with the corresponding eigenvalues of "H (i,i + 1). 
Then we confirmed that J c s/5\ and K e g/5 A are almost 
independent of 5a typically for 5a < 0.002. The constant 
values of J e g/5 A and K c g/5 A are shown in Fig. 4 for three 
phases (n=0, 1, and 2), which will be explained below. 

Because the effective coupling constants satisfy < 
K c ff/J c tf < 1, the ground state is the Haldane state for 
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Kaff/Jff 



Kaff/Jff 



Koff/Jff 



A, c (2,3) 



X c (l,l) 



Fig. 4. Effective bilinear interaction (J e ff) and biquadratic inter- 
action (K c ff) between spin clusters for small 5,4 in HDCO, HDC1, 
and HDC2 phases from top to bottom. The ratio K B g/J e g is also 
shown. 



small Sa- 3JS> In the Haldane state, each spin-1 degree of 
freedom is carried by a cluster-n rather than by a sin- 
gle spin. We call the state the Haldane DCn (HDCn) 
state. In the HDCn state with n > 1, the (n + l)-fold 
translational symmetry is spontaneously broken unlike 
the conventional Haldane state without STSB. Both the 
HDCO state for A > A c (0, 1) and the HDCoo state for 
A < A c (3,oo) are the Haldane states without STSB. In 
particular, the HDCoo state continues from the Haldane 
state (DCoo state) of the undistorted MDC mentioned 
in §2. 12 ) 



(a) 



(b) 



(e) 




Uniform Haldane 
X>\ c (t ,0) 



HDC1 



(c) fc 



(d) 



HDC2 



HDC3 



Uniform Haldane 
\<\ c (3,o) or 5 A » 1 



Fig. 5. Valence bond structures of the ground states of all phases 
for the MDC with type A distortion. A small filled circle repre- 
sents a spin with a magnitude of 1/2. An original spin with a 
magnitude of 1 is represented by two decomposed 1/2 spins in 
an open circle indicating the symmetrization. A valence bond is 
represented by a dashed oval. 



3.2 Connection to the strong distortion regime 

In the strong distortion regime of <5a —► 1 and small 
A, the three spins rfj-p Si, and t\ ' form a singlet clus- 
ter. Hence, the ground state is a state with spin gap and 
without STSB. This nature is common to the HDCO and 
HDCoo phases in §3.1. Furthermore, the HDCoo state 
is transformed into the HDCO state only by rearranging 
two valence bonds within each diamond unit, as shown 
in Figs. 5(a) and 5(e). Therefore, the strong distortion, 
HDCO, and HDCoo regimes are considered to be differ- 
ent parts of a single phase. The continuity of the three 
regimes will be confirmed by the numerical analysis dis- 
cussed in §3.3. In what follows, we call this phase the 
uniform Haldane (UH) phase as a whole. 

3. 3 Numerical phase diagram 

Under the periodic boundary condition, even in the 
parameter region where STSB takes place, the ground 
state of a finite chain is a superposition of the symmetry- 
broken states, and the translational symmetry is re- 
covered. Under the open boundary condition, however, 
one of the symmetry broken states is selected by the 
boundary effect. Therefore, we employ the DMRG cal- 
culation with the open boundary condition to determine 
the phase diagram for finite 6a- The DMRG calculation 
is carried out using the finite-size algorithm up to 288 
sites keeping 200 states in subsystems. We calculate the 
ground-state expectation values (T l ) and define the ef- 
fective spin magnitude TJ on the Z-th diagonal bond by 
Tt(T t + 1) = (T 2 ). A typical / dependence of T t is shown 
in Fig. 6 in each phase. With the increase in <5a, the trans- 
lational symmetry is recovered as expected. For finite 6a, 
the ground-state phase is identified from the periodicity 
in the oscillation of Tj. In the HDCn phase, the values of 
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Fig. 6. Profiles of T; for (a) A = 2.8 and (b) A = 2.62. 
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Fig. 7. System size dependences of AT at (a) A = 2.85 and (b) 
A = 2.62. The data are plotted against N~P'" where /3 and v 
are the critical exponents of the order parameter and correlation 
length, respectively, for the 2-dimensional (a) Ising and (b) 3- 
clock model. 



Ti follow the sequence 

...T S T L ---T L T S ,T L ---T L 



(T L > T S ). (7) 



Thus, we define the order parameter of the HDCn phase 
by AT = T L - T s . In DMRG, AT is measured at the 
sites closest to the center of the chain. 

The valence bond structures for the HDCn phases as 
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Fig. 8. Phase diagram of the MDC with type A distortion. The 
triangles indicate the position of the phase boundary for <5a = 0. 



well as the UH phase are shown in Fig. 5. We see the 
translational invariancc of period n + 1 in the HDCn 
ground state in contrast to the period- 1 invariance in 
the UH ground state. Hence, the Z n+ i STSB takes place 
at the HDCn-UH phase boundary. We expect that this 
transition belongs to the 2-dimensional (n + l)-clock 
model universality class. The system size dependence 
of AT for A = 2.85 is shown in Fig. 7(a) around the 
HDC1-UH phase boundary. Here, the data are plotted 
against N~ l3 / V with the order parameter critical expo- 
nent j3 = 1/8 and the correlation length critical expo- 
nent v = 1 for the two-dimensional Ising universality 
class. This shows that the critical value of <5a lies between 
0.0304 and 0.0308. A similar plot is shown in Fig. 7(b) 
for A = 2.62 around the HDC2-UH phase boundary, as- 
suming the critical exponents of two-dimensional 3-clock 
model (equivalently 3-state Potts model 38 )) with (3 = 1/9 
and v = 5/6. This shows that the critical value of 5a 
lies between 0.00822 and 0.00826. The critical points at 
other values of A are determined similarly. The results 
are shown in the phase diagram of Fig. 8. The error bars 
are within the size of the symbols. 

To confirm the consistency of the universality class, 
the finite size scaling plot for the order parameter AT is 
carried out. According to the scaling hypothesis, the 5a 
dependence of the order parameter AT of the finite size 
systems near the critical point should obey the finite size 
scaling law 39 ) 



ATN^ = f(N(6 A -6 c A r), 



(8) 



in terms of the scaled variables ATN^/ U and N{6a~ 5 a ) 1 ' 
and the scaling function f(x). In Figs. 9(a) and 9(b), 
ATA' 3 /" is plotted against N(6 A - 6\) v around the 
HDC1-UH and HDC2-UH phase boundaries assuming 
the Ising and 3-clock universality classes, respectively. 
The critical points 6 A = 0.0307 (Fig. 9(a)) and 0.008248 
(Fig. 9(b)) are chosen so that all data fall on a single 
universal scaling curve as well as possible. These plots 
are consistent with the expected universality class. 
The critical behavior at the HDC1-UH transition in 
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Fig. 9. Finite-size scaling plot of AT around the critical points, 
(a) Plot around the HDC1-UH phase boundary at A = 2.85. 
The Ising critical exponents v = 1 and /3 = 1/8 are assumed. 
The critical point is set at <5^ = 0.0307. (b) Plot around the 
HDC2-UH phase boundary at A = 2.62. The 3-state Potts critical 
exponents v = 5/6 and /3 = 1/9 38 ' are assumed. The critical 



point is set at 8% = 0.008248. 
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Fig. 10. Valence bond structures of the ground states of spin-1 
bilinear-biquadratic chain in the (a) Haldane phase and (b) dimer 
phase. The spins with a magnitude of unity represented by open 
circles are decomposed into two spin-1/2 degrees of freedom rep- 
resented by small filled circles. The valence bonds are represented 
by dashed ovals. The spins belonging to disconnected clusters in 
the dimer phase are connected by the valence bonds in the Hal- 
dane phase. 



our model should be compared with that of the S — 
1 bilinear-biquadratic chain at the Takhtajan-Babujian 
point. 40 ' 41 ) Both transitions are accompanied by Z%- 
STSB which contributes to the conformal charge by 1/2. 
For the HDC1-UH transition in our model, the rearrange- 
ment of valence bonds take place only within each dia- 
mond unit, as shown in Fig. 5(a) and 5(b). In contrast, 
in the S = 1 bilinear-biquadratic chain, the spins be- 
longing to disconnected clusters in the dimer phase are 
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Fig. 11. Phase diagram of the MDC with type B distortion with 
3N = 18. The triangles indicate the position of the phase bound- 
aries Ac (n,n+ 1) for <5b = 0. 12 ) 



connected by the valence bonds in the Haldane phase, as 
shown in Fig. 10. Apart from Z2-STSB, this is similar 
to the Gaussian criticality of the Haldanc-dimer transi- 
tion in the spin-1 alternating bond Hciscnberg chain that 
contribute to the conformal charge by l, 42-44 ) Therefore, 
the S = 1 bilinear-biquadratic chain at the Takhtajan- 
Babujian point is described by the conformal field theory 
with c = 1/2 + 1 = 3/2, while the HDC1-UH transition 
in our model is described by the c = 1/2 Ising conformal 
field theory. 

4. Ground-State Properties of the MDC with 
Type B Distortion 

In the case of type B distortion, the effective inter- 
action between the spins of two clustcr-n's separated by 
the dimer consisting of r, and t\ is ferromagnetic, be- 
cause both Si and Si+\ tend to be antiparallel to t\ '. 
Therefore, we expect the ferrimagnetic ground state with 
spontaneous magnetization quantized as m = l/(n + 1) 
per unit cell for small <5b in the range A c (n, n + 1) < 
A < X c (n— 1, n). We call this phase a ferrimagnetic DCn 
phase (FDCn phase). In contrast, the ground state for 
A < A c (3,oo) will remain in the Haldane phase, since a 
nonmagnetic gapped phase is generally robust against a 
weak distortion. For finite <5b , we determined the ground- 
state phase diagram by the numerical diagonalization for 
the system size 3iV = 18, as shown in Fig. 11. Among 
system sizes tractable by numerical diagonalization, only 
this size of 3A = 18 is compatible with all the ground- 
state structures with n = 0,1, and 2. As expected, the 
FDCn phases with m = l/(n + 1) arc found for these 
values of n. 

By inspecting numerical data for the 3N = 18 system, 
we also find other narrow ferrimagnetic phases between 
the FDCn and FDC(n + 1) phases with n = 0, 1, and 
2, although they are too narrow to be shown in Fig. 11. 
In order to investigate these phases in detail, we employ 
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(a) S B =0.2 

3N=18(ED) 
3N=72(DMRG)oo 



r 




N, is related to the pseudospins erf as 



2.6 



2.8 




0.5^ 



Fig. 12. Spontaneous magnetization for (a) <5b = 0.2 and (b) 
5b = 0.6. The triangles on the vertical axes indicate the values 
of the spontaneous magnetization m = l/(n + 1) in the FDCn 
phases. 



the DMRG calculation for 3iV = 72 keeping 120 states in 
each subsystem. Typical examples of the A dependence 
of spontaneous magnetization are shown in Fig. 12(a) 
for 5 B = 0.2 and Fig. 12(b) for 5 B = 0.6. Between the 
FDCn and FDC(n + 1) phases with n = 0, 1, 2, we find 
the partial ferrimagnetic phase in which the spontaneous 
magnetization varies continuously with A. The ferrimag- 
netic phase of this kind has been found in various frus- 
trated one-dimensional quantum spin systems. 45 ~ 51 ) In 
contrast, between the nonmagnetic phase and the FDC3 
phase, we find no partial ferrimagnetic phase for small 
5 B . 

This can be understood as follows: At A = X c (n, n+ 1), 
the cluster-n and cluster- (n + 1) can coexist. As stated 
above, it is physically evident that the effective magnetic 
interaction between the clusters is ferromagnetic. There- 
fore, we can restrict the states of each cluster to the max- 
imally polarized ground state with Sf = 1. Hence, the 
ground state of the whole chain is described by specifying 
the arrangement of cluster-n's and cluster-(n + l)'s. We 
map the two possible values of the length of i-th cluster, 
m = n and n, = n + 1, to two possible values of the 
spin-1/2 pseudospin, erf = 1/2 and erf = —1/2, respec- 
tively. Then, the total magnetization M is equal to the 
number of clusters N c . The total number of unit cells, 



N 



E 



N r 



Ay 

E< 

£=1 



(9) 

Therefore, the ground-state magnetization per unit cell 
m is given by 

Nr 1 

(10) 



(N) 



n + ^-a 



,jV c 



with a = J2i=i (°f ) /N c - The bracket (• • •} represents the 
ground-state expectation value. In the presence of 5b, 
the length of neighboring clusters can exchange through 
a second order process in 5b- This corresponds to the 
spin exchange in terms of pseudospins. In this case, the 
interaction between the pseudospins is approximated by 
the spin-1/2 XXZ Hamiltonian 



Hxxz 



E 



%xxz(«, * + !)> 



(11) 



HXXZ(M + 1) = Jf <Tl°! + l + Jf(« +1 + ciVf +1 ) 

(12) 

up to the second order in 5b ■ Here, further neighbor inter- 
actions are neglected. We estimate the effective exchange 
constants by comparing the energy spectrum of the pair 
Hamiltonian TLxxz(i, i + 1) with that of the correspond- 
ing pair of clusters as follows: 

(i) A = A C (0,1) 

Jf ~ -0.0395b 2 , 



Jf ~ 0.0875b • 



(ii) A = A C (1,2) 



(hi) A -A c (2, 3) 



TCff 



TCff 



Jf 



J 



Off 



-0.0082fe 
0.069<5 B 2 - 



-0.00295b' 

0.0185b 2 . 



(13) 



(14) 



(15) 



The details of the calculation are explained in Appendix. 
In all cases (i) ~ (hi), we hnd that the effective cou- 



pling constants satisfy the inequality 



|jf| 



< J 



elf 



< 



\Jf\- As is well known, the ground state of the spin- 
1/2 XXZ chain in this parameter regime is nonmagnetic 
and gapless in the absence of a magnetic field. Roughly 
speaking, AA = A — A c (n, n + 1) corresponds to the effec- 
tive magnetic field h e s conjugate to the total pseudospin 
2i°f) because the increase in A favors cluster-n over 
cluster-(n + 1); however, this correspondence should not 
be taken literally. A more precise argument is also given 
in Appendix. When AA takes a large negative value, 
the pseudospins are fully polarized downward to give 
(erf) = —1/2. This state corresponds to the FDC(n + 1) 
state with m = l/(n + 2). When h B s reaches the critical 
value h c i = — (\JJ_ \ + Jf), the magnetization starts to 
increase continuously until all spins are fully polarized 
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upward at the critical effective field h C 2 = \J]^\ + J% ■ 
This corresponds to the FDCn state with rn = l/(n+l). 
On the other hand, the magnetization jumps from to 
1/4 at the phase boundary between the Haldane phase 
and the FDC3 phase for small 8b- At this phase bound- 
ary, no finite size clusters coexist with cluster-3. There- 
fore, no pseudospin degrees of freedom can be defined. 
Consequently, no partial ferrimagnctic phase can be re- 
alized. In contrast, for larger values of 5b, we numerically 
find a partial ferrimagnctic phase between the FDC3 and 
UH phases. This would be ascribed to the contribution of 
other finite-length clusters with low lying energies which 
come into play through higher-order processes in 5b- 

5. Summary and Discussion 

We introduced two types of distortion, type A and type 
B, into the MDC with spins 1 and 1/2, and investigated 
the ground-state phases. The phase diagrams are charac- 
teristic of the type A and type B distortions, respectively. 
For the type A distortion, the effective interaction be- 
tween the cluster spins is antifcrromagnctic with bilinear 
and biquadratic terms. The numerically estimated val- 
ues of the effective couplings show that the DCn ground 
states are transformed into the HDCn ground states. The 
order parameters characterizing the HDCn phases are 
defined and the UH-HDCn phase boundaries are deter- 
mined using the DMRG data. From the valence bond 
structure of each phase, we expect that the UH-HDCn 
phase transition belongs to the universality class of the 
2-dimcnsional (n + l)-clock model. The finite size scaling 
plot of the order parameter is consistent with this iden- 
tification. For the type B distortion, the effective inter- 
action between the cluster spins is ferromagnetic. In ad- 
dition to the FDCn phases with quantized spontaneous 
magnetization m = l/(n, + 1), the partial ferrimagnctic 
phases are also found numerically between the FDCn and 
FDC(n+ 1) phases. A physical interpretation of the par- 
tial fcrrimagnetic phase is given for small 5b by mapping 
onto an effective pseudospin- 1/2 XXZ chain. 

Generally, the introduction of lattice distortion into 
a physical model increases the possibility that a corre- 
sponding material is realized. In the MDC, there are 
three types of distortion modes affecting the exchange 
interactions. Among them, the two types investigated in 
the present paper are of generic nature, because the lo- 
cal conservation laws that hold in the undistorted MDC 
are broken. This suggests that the observation of the ex- 
otic phenomena predicted in the present paper is pos- 
sible even if the corresponding material is not exactly 
described by the model Hamiltonians (1) and (2). 

If a distorted MDC material is synthesized, the dis- 
tortion may be controlled by, e.g., applying pressure. If 
the distortion is of type A, the Curie constant vanishes 
as the DCn ground state turns into one of the HDCn 
ground states. The magnetic susceptibility and magnetic 
specific heat will have an activation-type temperature 
dependence with activation energy proportional to the 
effective coupling between the cluster spins, which is of 
the order of 5 a 2 - These HDCn phases are not realized 
if the distortion 5a exceeds ~ 0.03 even in the most ro- 
bust case of n = 1. In a real material, the STSB in the 



valence bond structure manifests itself as a magnetic su- 
perstructure. It is also possible that it is accompanied 
by a lattice superstructure of corresponding periodicity 
if the spin-lattice coupling is present. Therefore careful 
measurements of magnetic and lattice superstructures 
would help with the observation of HDCn phases with 
1 < n < 3. 

On the other hand, if the distortion of the material 
is of type B, the ground sate is ferrimagnetic. At low 
but finite temperatures, however, the spontaneous mag- 
netization vanishes owing to the one-dimensionality As 
a precursor of ordering at T = 0, the low-temperature 
magnetic susceptibility should diverge as T~ 2 with a co- 
efficient proportional to the effective couplings 5b be- 
tween the cluster spins. 52 ~ 54 ) This means that even a 
weak magnetic field of the order of H ~ T 2 /5b 2 derives 
the finite-temperature magnetization up to the value of 
the ground-state spontaneous magnetization. This en- 
ables the experimental estimation of the spontaneous 
magnetization in real materials. The quantized ferrimag- 
netic behavior should be observed for wide ranges of the 
parameters A and 5b as shown in Fig. 11, and should 
be easily observed if an appropriate material is synthe- 
sized. The partial ferrimagnctic phases are limited to 
narrow intervals of the parameters 5b and A. Therefore, 
these can only be observed as a temperature- independent 
crossover between two quantized ferrimagnctic behaviors 
with careful exclusion of the thermal effect. 

We have demonstrated that various exotic ground 
states and phase transitions between them are realized 
in the distorted MDC with spins 1 and 1/2, which has 
a strong frustration. The physical pictures of these phe- 
nomena have become clear. This is made possible because 
the ground state of the undistorted MDC is known ex- 
actly. Therefore, we expect that our model may provide 
a means of understanding the similar exotic phenomena 
realized owing to the interplay of spin ordering, quantum 
fluctuation, and strong frustration in more general frus- 
trated quantum chains on a firm ground. For example, 
partial fcrrimagnetic phases are found in various one- 
dimensional frustrated quantum spin models. 45 ~ 51 ) How- 
ever, some of them are only numerically confirmed and 
no physical explanation has been given so far. We hope 
that the present study paves the way to the general un- 
derstanding of these partial ferrimagnctic states. 
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Appendix 

The Hamiltonian Hb with the type B distortion is 



rewritten as 



Hb — T~Ln 



sn, 



where 



N 



5H = 5 B J2{Si + S l+1 ){T< i 1) - 



.(2)1 



(A-l) 



(A-2) 



1=1 



For small 5b, the ground state around A = A c (n, n + 1) 
consists almost entirely of cluster-n's and cluster- (n + 
l)'s. Hence, as a good approximation, we consider Kb 
only in the restricted Hilbert space where each state 
involves no clusters except for clustcr-n's and cluster- 
(n + l)'s. Under the fixed cluster number 7V C in this 
Hilbert space, "Ho is equivalent to the following effective 
Hamiltonian expressed in terms of pscudospin operators: 

H eS = E° G (n+l;\)Y,(--a* 



+ E° G (n;X) 



E 



2 + ° l 



(A-3) 



where of = 1/2 and of = —1/2 correspond to n% = n 
and rii = n + 1, respectively. E G (n;X) is the ground- 
state energy of a cluster-n and a dimer in the absence of 
distortion, and is given by 

^ n 3A 

T - T' (A4) 



E Q {n;X) = E{2n + l) + 



where E(2n + 1) is the ground-state energy of the spin-1 
antifcrromagnctic Heisenbcrg chain with length 2n+l. 12 ' 
The application of 5H to the unperturbed ground state 
transforms one of the TJ = bonds to a 1} = 1 bond or 
vice versa. Then the resulting states contain clusters with 
lengths less than n or greater than 2n. Since these states 
are outside the restricted Hilbert space, no correction to 
the ground-state energy is present within the first order 
in 5b • Hence, the lowest-order correction is of the order 
of 8b 2 - Up to the second order in S-q, the effective pseu- 
dospin Hamiltonian is given by 



H e & = E G (n + 1; A, S B ) ^ ( \- <*■ 



t=i 



Ec(n; A,<5b) 



Ay 

E 

i=l 



u 



XXZ; 



(A5) 



where Eq(ji; A, 6b) is the ground-state energy of a 
cluster-n and a dimer including the second order cor- 
rection in 5b ■ This is simply expressed as 



N c 



H 



cir 



N C E G + AE G J2 of + n XX z, 



(A-6) 



with 



Eg = ^(E G (n + 1; A, 5 B ) + E G (n; A, 5 B )), (A-7) 



AS G = E G {n; A, <5 B ) - E G {n + 1; A, 5 B ) 



(A- 



The effective coupling constants J° and J^_ in Hxxz 
are also of the second order in 5b- We determine J° and 
J^_ so as to reproduce the low-lying energy spectrum of a 
pair of cluster-n's by that of two-pseudospin Hamiltonian 

H cS {i, i + l) = 2E G + AE G {a* + of +1 ) + Hxxz(i, i + 1) 

(A9) 



In each of the subspaces of + a 



i+i 



±1, of 



cf+i = C = ±1/2. Therefore, the Hilbert space is one- 
dimensional and the eigenvalue of H e s(i,i + I) is sim- 
ply E aa = 2E G + 2AE G <7 + Jf/A with a = ±1/2. 
In the subspace of + o. 



i+i 



0, the Hilbert space is 
two-dimensional and the eigenvalues of 'Hoff(*,* + 1) are 
E± = 2E G - Jf/A ± Jf/2. 

The original Hamiltonian of the cluster consisting of a 
cluster-n and a cluster-n' is the distorted diamond chain 
with length n + n' . 

n+n' + 1 

H(n + n')= J2 [(1 + <5b)^t / (1) + (1 + ,5b)t| 1) ^ +1 
i=i 

+ (I - 5 b )S 1 t\ 2) + (1 - 5 b )t\ 2) S 1+1 + At, (1) t, (2) " , 

(A-10) 

We denote the a-th eigenvalue of ~K{n + n 1 ) as E(n + 
n';a). Comparing the corresponding expression for the 
eigenvalues, we find 

T0ff 

E(2n; 0) = Ex i = 2E G + AE G + -*- (A-ll) 



£;(2n + 2;0) = E_ 



E(2n+l;0) = E_ =2E Q 



Teff 

2£ G -A£ G + ^- (A12) 



7-cff Teff 

J z J + 



(A-13) 



reff Teff 

E(2n + l;l)=E+ = 2EG--%- + ^-. (A-14) 

Solving these sets of equations, with respect to J| ff and 

(A-15) 



J^, we find 



Jf = E(2n + 1; 1) - E(2n + 1; 0), 

Jf = 2[E(2n + 2; 0) + E(2n; 0) 

-E(2n+l;l)-E{2n+l;0)}. (A-16) 

Note that the rhs's of (A-15) and (A-16) vanish for 
5 B = 0. We numerically evaluated E(2n; 0), E(2n+1; 0), 
£(2n + 1; 1), and E(2n + 2; 0) at A = A c (n, n + 1) (n = 
0, 1, 2) for small 5b- Using these values in eqs. (A-15) and 
(A-16), we determined JJ_ and J° as eqs. (13)-(15). 

For the whole MDC, the ground-state energy is written 
as 



E = N C E G + N c AE G a + N c e xxz (a) 



(A-17) 



where exxz(c) is the ground-state energy per site of a 
magnetized spin-1/2 XXZ chain with (of) = a. The 
number of unit cells, TV, of the original MDC is given by 
the expectation value of eq. (9) as N = N c (n + | — a) . 
Therefore, we have 



E n 



N 



[E G + AE c a + exxzW) - (A-18) 
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Fig. AT. Relationship between c and AA for n = 0, 1, and 2. 



Minimizing this with respect to a with fixed N, we find 

AX=(n+ 3 --^ ^^1 + exxzta), (A-19) 

where AA = A — A c (?i, n + 1; 5b) and A c (n, n + 1; 5b) is 
defined by 

(n + 2)£ G (n; A c , fa) - (n + l).Efa(n + 1; A c , fa) = 0. 

(A-20) 

To simplify the calculation, we replace exxz(c) by the 
ground-state energy of the spin- 1/2 XY chain exy = 
— (J^/ir) cos7rcr, because \J^\ is substantially larger 
than I JJ ff | in all cases. Then we find 



AA 

1W 



3 A 1 

n . -) cr J sin ira cos ira. (A-21) 

2 / 7T 



This relation is plotted in Fig. A-l for n = 0, 1 and 2. It 
is clear that a continuously increases from —1/2 to 1/2 
with an increase in A. 
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